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We theoretically study the carrier-envelope phase dependent inversion generated in a two-level system by 
excitation with a few-cycle pulse. Based on the invariance of the inversion under time reversal of the exciting 
field, parameters are introduced to characterize the phase sensitivity of the induced inversion. Linear and 
nonlinear phase effects are numerically studied for rectangular and sinc-shaped pulses. Furthermore, analyt- 
ical results are obtained in the limits of weak fields as well as strong dephasing, and by nearly degenerate 
perturbation theory for sinusoidal excitation. The results show that the phase sensitive inversion in the ideal 
two-level system is a promising route for constructing carrier-envelope phase detectors. 
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1. INTRODUCTION 

The two-level system is a fundamental and widely 
used model to describe the interaction of electromagnetic 
waves with matter. The temporal evolution of the sys- 
tem is usually treated within the framework of the rotat- 
ing wave approximation (RWA), where the driving elec- 
tric field enters the equations of motion only through its 
complex envelope and center frequency^ For a complete 
description of the electric field, also the carrier-envelope 
(CE) phase has to be taken into account, specifying the 
position of the envelope with respect to the rapidly oscil- 
lating carrier wave. The RWA breaks down for strongly 
driven systems,— giving rise to new effects which not only 
depend on the pulse envelope and carrier frequency, but 
also on the CE phase. 

The phase sensitive dynamics of the driven two-level 
system beyond the RWA has been the topic of several 
papers. The phase dependence of the inversion was care- 
fully examined for a sinusoidal excitation, serving as a 
model for the interaction of atoms and molecules with 
continuous wave laser radiation^ - — The discussion was 
extended to rectangular pulses, obtained by abruptly 
switching on and off the sinusoidal excitation. The phase 
sensitivity of the inversion was investigated and exper- 
imentally demonstrated in the radio frequency regime 
by exciting the anticrossing of the potassium 21s - 19f 
states^ 

Following the arrival of laser pulses consisting of only 
a few optical cycles, there has been considerable in- 
terest in phase sensitive effects in the pulsed optical 
regime^, and various approaches have been used for de- 
tecting the CE phase/frequency»i2.~— In this context, the 
CE phase dependent emission of two-level systems^! and 
semiconductors— interacting with few-cycle pulses has 
been theoretically investigated, and the effect has exper- 
imentally been observed in GaAs.- — So far, less attention 
has been given to the CE phase dependence of the inver- 
sion, with a few exceptions studying the interaction with 



Gaussian pulsesi 20 ' 21 

While for sinusoidal excitation, the generated inversion 
shows a CE phase dependence even for weak fields, the 
phase sensitive dynamics relies completely on nonlinear 
effects for pulsed optical excitation. In this paper, both 
linear and higher order phase dependent inversion effects 
in two-level systems are theoretically investigated, us- 
ing analytical approximations and numerical simulations 
with properly chosen test pulses. In addition, the in- 
fluence of dephasing on the phase sensitivity is studied. 
General properties of the steady state inversion are dis- 
cussed, and approximate expressions are derived in the 
linear and the nonlinear regime. The paper is organized 
as follows: In Section 2, the equations of motion for an 
excited two-level system are given in a fixed and a ro- 
tating reference frame. In Section 3, general properties 
of the steady state inversion are discussed. In Section 
4, the phase dependent inversion in the weak field limit 
is analytically examined, and in Section 5, the discus- 
sion is extended beyond the linear regime for rectangular 
and sinc-shaped pulses, using numerical simulations and 
nearly degenerate perturbation theory. In Section 6, the 
influence of the phase relaxation is studied based on the 
strong dephasing approximation and numerical results. 
We conclude in Section 7. 



2. EQUATIONS OF MOTION 

A two-level system is characterized by its dipole ma- 
trix element d and resonance frequency uiba — ^fba 

/.'; — E a )/h, where E a and Eb are the eigenenergies as- 
sociated with the low and high energy states \a) and \b), 
respectively. Dissipative effects, which arise due to the 
interaction of the ideal two-level system with its envi- 
ronment, can be taken into account in a statistical ap- 
proach. The density matrix is represented by the com- 
ponents of the Bloch vector s. The components si and 
S2 are related to the real and imaginary parts of the off- 
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diagonal density matrix elements by p a b = (si + is 2 ) /2, 
and — p aa = s 3 = w is the population inversion. 
In this paper, the relaxation processes are described by 
phenomenological parameters, the energy relaxation rate 
71 = 1 / T\ and the dephasing rate 72 = I/T2 . Frequently, 
relaxation is dominated by processes which lead to a de- 
struction of the phase coherence in the quantum system 
without affecting the inversion, resulting in a dephasing 
time T2 which is much shorter than the energy relaxation 
time T\. For example, this is typically the case in a gas 
due to collision broadening. 22 Thus, in the following, we 
set 71 = 0, assuming that the energy relaxation processes 
are slow compared to the interaction time with the field, 
whereas we do allow for dephasing processes occurring on 
a time scale comparable to the duration of the exciting 
pulse. 

Assuming linear polarization of the exciting field and 
a vanishing static dipole moment, the dynamics of the 
system is described by the Bloch equations, 1 



S\ = -U ba S 2 - 72 Si, 

h = Wfe a Si + 2fLs 3 - 72S2, 

S3 = -2fls 2 . 



(1) 



The overdot denotes a time derivative. The electric field 
E (t) is parametrized in terms of the instantaneous Rabi 
frequency f2 (t) = dE (t) / h, which can be written as 



O (t) = Q R [e (t) exp (-kj c t + i^cE) 
+e* (t) exp (iui c t - h^ce) ] A 



(2) 



where the asterisk denotes the complex conjugate. Here, 
e (t) is the normalized, in general complex envelope func- 
tion, ui c = 27r/ c is the carrier frequency, and (/>ce is the 
CE phase. In this paper, we refer to Or = dEo/h as the 
(peak) Rabi frequency, with the maximum value of the 
electric field envelope Eq. 

For some applications, it is convenient to transform the 
Bloch equations, Eq. ([T]), into a rotating reference frame. 
The transformation between the Bloch vector s and the 
vector in the rotating frame u = [u, v, w) , where T de- 
notes the transposed vector, is given by u (t) = A (t) s (t), 



cos (ui ba t) sin (u)b a i) 
A (t) = I - sin (uj ba t) cos (uibat) 




(3) 



For 71 = 0, we obtain the equations of motion in the 
rotating frame 



u = 20 sin (ui ba t) w — 72W, 

V = 2f2 COS (uibat) w — 72V, 

= —20 sin (ujbat) u ~ 257 cos (uibat) v. 



w 



(4a) 
(4b) 
(4c) 



3. STEADY STATE INVERSION 



We consider a two-level system excited by an electric 
field pulse, parametrized in terms of f2(£), and neglect 



energy relaxation, i.e., 71 = 0. During interaction with 
the field, the inversion in the system evolves and reaches 
a steady state value w s , staying constant after interaction 
with the pulse. In this paper, we investigate the CE phase 
sensitivity of the steady state inversion, which is readily 
accessible to experiments. We assume that the dipole 
moment of the system is initially vanishing, correspond- 
ing to the initial condition si = s 2 = and S3 = wq. This 
is always fulfilled if the system starts from equilibrium, 
where wq corresponds to the equilibrium inversion. 



A. Inversion invariance under time reversal of the field 

Strong excitation of the two-level system results in the 
breakdown of the RWA, and of related features like the 
area theorem* 2 - By way of contrast, the property dis- 
cussed in the following is even valid in the strong field 
limit. 
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FIG. 1. Dynamics of the two-level system for excitation with 
a field (t) (black curve) and the corresponding time-reversed 
field Q (— t) (gray curve), (a) Time dependent electric field, 
parametrized in terms of the normalized Rabi frequency, (b) 
Evolution of the inversion for a Rabi frequency Or = 1.3cj c , 
a transition frequency u)i a — 1.7 u> c , and a dephasing rate 
72 = /c/10. 

We take the Bloch equations, Eq. (f?]), assuming 71 = 
and an initially vanishing dipole moment. Under these 
conditions, the steady state inversion w s after interaction 
with an arbitrary pulse f2 (t) is the same as for excita- 
tion with the time-reversed pulse fl (— t). The proof is 
given in Appendix A. This feature is remarkable, since 
the temporal evolution of the system from its initial state 
is completely different in both cases, and should not be 
confused with the time-reversed dynamics of the system. 
In Fig. [TJ the invariance of w s with respect to time re- 
versal of the field is illustrated by means of an example. 
The time dependent field of an exciting pulse is shown 
in Fig. [T]Ja), as well as the corresponding time- reversed 
field. The pulse shape is in principle arbitrary; the only 
condition is that the pulse has a finite energy, and thus 
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FIG. 2. Bloch vector trajectories for excitation of a two-level 
system with a driving field (black curve) and the time-reversed 
field (gray curve). The fields and system parameters are the 
same as for Fig. [1] 



the field approaches zero for t —¥ ±oc. For a propa- 
gating electric field, we furthermore require a vanishing 
dc component. This condition, however, is not crucial 
for the inversion invariance discussed here. In Fig. [TJb), 
the time dependent inversion is displayed, starting at the 
initial value wo = — 1. The evolution of the inversion 
is completely different for excitation with the field and 
the time-reversed field, but still, the same steady state 
value is reached in both cases. This can also be seen in 
the Bloch sphere representation. In Fig. [2j the trajecto- 
ries of the Bloch components u, v and w in the rotating 
frame are shown, with the south pole of the sphere cor- 
responding to the initial state. The Bloch vectors evolve 
in completely different ways, but the same steady state 
inversion is reached. In the presence of dephasing, the 
case which is shown here, the trajectories approach the 
same point for t — > oo, since the dipole moment decays, 
i.e., u and v tend to zero for t — > oo. By way of contrast, 
for 72 = 0, u and v approach different steady state values 
for excitation with the field and the time-reversed field, 
respectively. Thus, for t — > oo, the two trajectories still 
reach the same value for w, but not for u and v. 



B. Steady state inversion for symmetric pulses 

In Sections 5 and 6, we discuss the two- level dynam- 
ics for excitation by sinc-shaped and rectangular pulses, 
which are chirp-free and have symmetric intensity en- 
velopes. For such pulses (or more generally for pulses 
with £ (i) = e* (—£)), the time reversal of the field cor- 
responds to a mere sign change of the CE phase, as can 
be seen from Eq. (5). As a consequence of the time re- 
versal invariance discussed above, we have w s (</>ce) = 
w s (— 0ce) for these pulses. Due to the centrosymme- 
try of the two-level system, the inversion is furthermore 
invariant with respect to a sign change of the electric 



field, corresponding to a t shift in the CE phase, i.e., 
w s (4>ce) — w a (<^ce + 7r) . Thus, the CE phase depen- 
dent occupation probability can be represented by a it 
periodic Fourier series. If w s (4>ce) = w s (— <^>ce) holds, 
only the even (i.e., cos) terms remain. For non-resonant 
excitation and non-excessive field strengths, the series 
can be truncated after the first order term, and we ob- 
tain 

w s (</>ce) = w s + A cos (2</> C e) ■ (5) 

The parameters used here are the CE phase averaged 
inversion, 

w s = [w s (0) + w s (tt/2)] /2, (6) 
and the modulation amplitude 

A = [to. (0) - w s {tt/2)} /2, (7) 
with -1 < w s , A < 1. 
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FIG. 3. Phase dependent steady state inversion w a after in- 
teraction with a two-cycle sine pulse, for Rabi frequencies 
f^R, = 1.5 Ldc (solid curve) and JIr, = 2.3u c (dashed curve). 
The transition frequency is w& a = 1.5 cj c , and the dephasing 
time is T 2 = 10// c . 

Fig. [3] shows the phase dependence of w s after inter- 
action of the two-level system with a sine pulse. For 
Or = 1.5 cj c , the phase dependence of w s can be well ap- 
proximated by Eq. (JS]). For f2pi = 2.3 u> c , higher order 
terms in the Fourier expansion result in a bump around 
the symmetry point. For excitation near the two-level 
resonances ujh a — (2n, + 1) u) c , n = 0, 1, 2, . . ., the approx- 
imation Eq. {I} already breaks down for smaller Rabi 
frequencies, i.e., more terms in the Fourier series are nec- 
essary for a full description. 

4. WEAK FIELD LIMIT 

As mentioned in Section 1, linear phase sensitive effects 
have been observed for custom tailored microwave and 
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radio frequency pulses, while for laser-generated optical 
pulses, a CE phase dependence can only be obtained in 
the nonlinear regime. The phase dependent dynamics of 
an excited two-level system can be treated analytically 
by a time dependent perturbation series expansion, with 
the first order approximation describing the evolution in 
the linear regime. We use the Bloch equations in the 
rotating reference frame, Eq. Q. As initial condition 
at ( -> — oo, we choose u = (0, 0, Wq) T . By formally 
integrating Eqs. (|4"a)l and (|4bl) and inserting the results 
into Eq. (I4cl) , we arrive at an integral equation for w (t) , 



w (t) = w Q 



dt' 



dr COS (uJbaT) 



The power spectrum |fi (uj)\ is only phase independent 
if e (uj) = for u) < —uj c , which, for chirp-free pulses, 
corresponds to the condition e (uj) — for \uj\ > uj c f^ 

Formally, Eq. (jHJ is a Volterra integral equation of 
the second kind, as can be shown by exchanging the or- 
der of integration. The weak field approximation, Eq. 
([9]), is thus equivalent to a first order Neumann series 
expansion^ yielding (w a — wq) oc fl 2 ^ and also A oc fl 2 ^ 
for pulses with phase dependent power spectra. Thus, 
the amount of generated inversion and its phase depen- 
dence scale linearly with the peak intensity of the pulse. 
This is different for pulses with phase independent power 
spectra, where in the weak field limit, the phase depen- 
dence of w s is provided by the second order term. In this 



x exp (-7 2 t) n (t r ) Q (t' -t)w (t' - t) . (8) case > we obtain A oc (or A oc ft| for two-level sys- 



In the weak field limit, the time dependent inversion 
changes only slightly with respect to its initial value. The 
first order approximation (which is second order in Or) is 
obtained by inserting the Oth order solution (t) = wq 
into the right-hand side of Eq. ([5]) , 



w {1) (t) = w 



/>oo 

1 - 4 / dr cos (uj ba r) exp (-72T) 
Jo 



dt'0,(t')Q(t' -t) 



(9) 



The result for the steady state inversion w s = w (t — » 00) 
is favorably expressed in terms of the Fourier transform 



Q(u)= dttt (t) exp (iwt) , 



yielding 



(i) 



Wq 



1 



duj \n(u>)\ 2 H(u) 



(10) 



(11) 



Thus, in the linear regime, the relationship between w s 
and the power spectrum of the pulse |fi (uj)\ 2 can be de- 
scribed in terms of a spectral filter function H (w), which 
is given by 



H(uj) 



72 (uj 2 + 7 | + UJ 2 a ) 



(12) 



For 72 <C uJba, we can make the approximation 

7T (U) - UJ ba ) + 7| 

and in the limit of no dephasing, 72 — > 0, we obtain 

H(u) = 25 (u> - u> ba ) . (14) 

The Fourier transform of fl (t), Eq. ((TUJ), can with Eq. 
(|2|) be written as 

Q (uj) = £1 R [e (uj - uj c ) exp (i^cE) 

+e*(-w-w c )exp(-i(/>cE)]/2. (15) 



terns with a static dipole moment^!), i.e., the generated 
inversion becomes phase insensitive for weak fields. 

In the microwave regime, pulses with the desired en- 
velopes and CE phases can be custom tailored. For ex- 
ample, pulses with durations far below one carrier cycle 
are available, giving rise to linear CE phase dependent 
processes^ For the experimental demonstration of the 
phase sensitive inversion, pulses in the radio frequency 
regime with approximately rectangular envelopes were 
applied, 8 which have a phase dependent power spectrum. 
For such pulses, linear filter effects can be used to deter- 
mine the CE phase, and a phase dependent inversion w s 
is obtained even in the weak field limit. Strictly speaking, 
this is also true for Gaussian or hyperbolic secant pulse 
shapes, often used as a model for laser pulsesi 20 i 21 ' 25 In 
contrast, for optical pulses emitted by current laser sys- 
tems, the power spectrum is phase independent. As a 
consequence, a CE phase dependent inversion w s can only 
be observed in the nonlinear regime. 



5. SYSTEMS WITHOUT DISSIPATION 

In the following, we investigate the phase dependence 
of the steady state inversion for interaction with rect- 
angular and sine pulses, serving as typical model pulses 
for excitation in the radio frequency and optical regime, 
respectively. In this section, we neglect dissipation, i.e., 
we set 71 = 72 = 0. The rectangular pulse corresponds 
to a sinusoidal field of finite duration and is thus closely 
related to continuous wave excitation, which plays an im- 
portant role in atomic and molecular physics and quan- 
tum chemistry^— This particularly basic pulse shape fea- 
tures a piecewise constant envelope. The periodicity of 
the field can be exploited for approximately solving the 
equations of motion, and basic features of the solutions 
can be derived. In the radio frequency regime, rectan- 
gular pulse shapes can be approximately generated by 
switching on and off a sinusoidal field, and they were 
used for the experimental demonstration of the phase 
sensitive inversion^ As already pointed out, rectangular 
pulses have a phase dependent power spectrum, yielding 
phase sensitive effects even in the linear regime. Such 
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pulses are currently not available in the optical regime 
from laser sources. Sine pulses, on the other hand, fea- 
ture a rectangular, phase independent power spectrum, 
and have proven to be a fairly good description of few- 
cycle laser pulses^ 

A. Excitation with rectangular pulses 
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FIG. 4. Rectangular two-cycle pulse: Power spectrum for CE 
phases (solid line) and it/2 (dashed line). 

For rectangular pulses, the envelope in Eq. ([2]) is given 
by e (t) = 1 for -T/2 < t < T/2, and e (*) = other- 
wise, where T is the pulse duration. In order to avoid a 
non-propagating dc component for arbitrary CE phases 
</>ce, the pulses must contain an integer number of op- 
tical cycles, i.e., T is an integer multiple of the carrier 
period T = l// c . For </> C e ^ (ra+ 1/2)tt, n E Z, the 
field strength exhibits discontinuities, giving rise to a CE 
phase dependent power spectrum. In Fig. [4J the power 
spectrum |0 (w)| of a rectangular two-cycle pulse is dis- 
played. The CE phase dependence of the power spectrum 
becomes also clear from the spectral properties of the 
pulse. Since the pulse spectrum does not vanish iden- 
tically for \uj\ > w c , the two components in Eq. (fT5j). 
centered around w c and — w c , have a region of spectral 
overlap, resulting in a CE phase dependent interference. 

In the following, we investigate the phase sensitivity 
of the inversion by numerically solving the Bloch equa- 
tions, Eq. (JT|). Fig. [5] shows the average inversion w s 
and the modulation amplitude A as a function of Or 
and ujba- In Fig. [Sfa) and (b), w s is displayed for excita- 
tion with rectangular single-cycle and two-cycle pulses, 
respectively. For weak fields, the inversion is largest for 
near-resonant excitation. This resonance is more pro- 
nounced for the two-cycle pulses, due to their narrower 
spectrum as compared to single-cycle pulses. For Rabi 
frequencies approaching or even exceeding w c , regions 
with strong inversion can also be found for off-resonant 
excitation due to higher order transitions. A is shown 



in Fig. [5jc) and (d), again for excitation with rectangu- 
lar single-cycle and two-cycle pulses. During interaction 
with the pulse, the inversion performs oscillations due to 
carrier-wave Rabi flopping, 2 and the steady state value 
after interaction with the pulse w s strongly depends on 
the pulse and system parameters. As a consequence, for 
w s as well as A, regions with positive and negative values 
alternate. For rectangular pulses, the generated inversion 
exhibits a phase dependence even in the linear regime, as 
discussed in Section 4. We introduce the relative modu- 
lation amplitude or modulation depth 

5 = A/ (w s - w ) , (16) 

which is helpful to discuss the phase sensitivity in the 
weak field limit. A value of 5 = ±1 indicates maximum 
possible CE phase modulation of the generated inversion, 
(5 = indicates no modulation. While the amount of 
generated inversion vanishes for low field strengths, the 
relative phase sensitivity of w s does not disappear in the 
case of phase dependent power spectra, yielding a finite 
value for 6 in the weak field limit. 

For interaction with a rectangular pulse shape, 6 does 
not depend on the pulse duration T. This is a conse- 
quence of the periodic excitation, as shown in Appendix 
B, and is only true in two-level systems without dissi- 
pation. An approximate solution beyond the weak field 
limit can be obtained by using almost degenerate pertur- 
bation theory^ - — where the perturbation parameter is 
the normalized Rabi frequency Or/£J c . The zeroth or- 
der perturbation theory, corresponding to the rotating 
wave approximation,— does not contain any CE phase 
dependence and thus fails in predicting phase sensitive 
effects— Following the procedure in Ref. 28, we obtain 
the first order perturbative result 

s = K„ - gg) K« + oj c ) + uj c ni 

P (w 2 fl + U; 2 ) (uj ba + CJ C ) - W c fi^ ' 

see Appendix C. Fig. [5] shows the perturbative approx- 
imation and the numerical solution for the modulation 
depth. The perturbative result is in good agreement 
with the exact solution, especially for moderate field 
strengths and moderate detuning, where the almost de- 
generate perturbative treatment converges best.— For 
small Rabi frequencies Or, w s approaches zero, but the 
relative phase sensitivity of the inversion does not disap- 
pear, i.e., the modulation depth stays finite. In the limit 
Or — > 0, where the perturbation theory becomes exact 
and coincides with the weak field approximation, Eqs. 
(ITU and flU}, we obtain 

5 v = {u 2 ba -ujl)/{ul a +ul) (18) 

(except for odd multiphoton resonances, where 8 becomes 
degenerate for Or — > and the approximation breaks 
down). Thus, in the weak field limit, S > for exci- 
tation below resonance, and 6 < for excitation above 
resonance. For resonant excitation (oj c = ujba), where 




R c R c 

FIG. 5. Average inversion w s and modulation amplitude A as a function of the Rabi frequency f2a and the transition frequency 
uita in units of uj c . Displayed is w s after interaction with rectangular (a) single-cycle and (b) two-cycle pulses, and A after 
interaction with rectangular (c) single-cycle and (d) two-cycle pulses. 
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FIG. 6. Modulation depth 5 as a function of the Rabi fre- 
quency Qr and the transition frequency uita in units of ui c : 
Perturbative (dashed curves) and exact numerical result (solid 
curves) . 



the population transfer reaches its maximum, we obtain 
(5 = 0, i.e., the inversion does not show a phase sensitiv- 



ity. 

As already mentioned, the RWA does not predict any 
phase dependence in the inversion. This might be coun- 
terintuitive, especially for pulses with CE phase depen- 
dent power spectra, where the generated inversion de- 
pends on the CE phase even in the weak field limit. A 
straightforward explanation is that the validity of the 
RWA is restricted to near-resonant excitation^ while 
on the other hand, the modulation depth in Eq. (|T5|) 
indicates a phase sensitivity of the inversion just for ex- 
citation off resonance. 



B. Excitation with sine-shaped pulses 

Sine pulses are characterized by the full width at 
half maximum (FWHM) value T of the intensity enve- 
lope and the Rabi frequency J7r . The normalized en- 
velope function e(t), introduced in Eq. @, is given by 
e(t) = sine (i/r) = sin (t/r) / (t/r). The FWHM pulse 
duration is T = 2.783r. These pulses have a rectangular 
spectrum, extending from f c — 1/ (2itt) to f c + 1/ (2ttt) . 
We require r > 1/ (2irf c ) in order to avoid unphysical 
dc components. Under this condition, the power spec- 
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FIG. 7. Average inversion w s and modulation amplitude A as a function of the Rabi frequency Qn and the transition frequency 
uita in units of co c . Displayed is w s after interaction with sinc-shaped (a) single-cycle and (b) two-cycle pulses, and A after 
interaction with sinc-shaped (c) single-cycle and (d) two-cycle pulses. 



trum is CE phase independent, see discussion after Eq. 
(fl"5|) . Thus the inversion does not exhibit a CE phase 
sensitivity in the linear regime. We note that Gaussian 
or hyperbolic secant pulse descriptions feature unphysi- 
cal dc components and a phase sensitive power spectrum. 
This can lead to erroneous results, especially in the weak 
field limit or for strong dephasing, see Section 6, where 
the remaining phase sensitivity is due to the phase de- 
pendent power spectrum, as opposed to nonlinear electric 
field effects. 

The Bloch equations, Eq. (pTJ) , are numerically solved. 
In Fig. the average inversion w s and the modulation 
amplitude A are displayed as a function of fipt and uiba- 
In Fig. Efa) and (b), w s is shown for excitation with 
sinc-shaped single-cycle and two-cycle pulses, and the 
corresponding A is displayed in Fig. [TJc) and (d). As 
in the case of excitation with rectangular pulses, see 
Fig. [5J regions with positive and negative values alter- 
nate for both w s and A. The inversion w s is largest 
for near-resonant excitation, and the resonance is more 
pronounced for two-cycle pulses, which have a narrower 
spectrum. Higher order transitions become relevant if 
JIr approaches or even exceeds w c , leading to regions 
with strong inversion also for off-resonant excitation. As 



discussed in Section 4, the steady state inversion shows 
no phase dependence in the linear regime, i.e., the modu- 
lation depth 5, Eq. (|16|) . approaches zero for sine pulses. 
A significant phase sensitivity can only be found for field 
strengths in the regime of carrier- wave Rabi flopping, 
where nonlinear effects play a significant role. For sinc- 
shaped pulses, which have smooth envelopes unlike rect- 
angular pulses, the phase sensitivity of w s depends con- 
siderably on the pulse duration, since the phase depen- 
dence of the field is more pronounced for shorter pulses. 
Thus, higher field strengths are necessary to obtain a sig- 
nificant modulation amplitude A for excitation by two- 
cycle pulses as compared to single-cycle pulses. 



6. INFLUENCE OF DEPHASING 

In the following, we investigate the influence of phase 
relaxation processes on the CE phase dependent inver- 
sion. In the limit of strong dephasing, an approximate 
expression can be derived. We start from the inte- 
gral equation Eq. ([5]). Since for strong dephasing the 
Rabi oscillations of the inversion are significantly damp- 
ened and the kernel decays rapidly, we can approximate 
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w (t' — t) w w {t'). By differentiation with respect to t, 
we obtain the first order differential equation 

[>00 

w{t) = -Aw(t) I dr cos(w ha r) exp(-7 2 r) fl(i') - r) 
Jo 

(19) 

with the solution 



w (t) = Wq exp 



poo 

-4/ dr cos(w ba r)exp(-7 2 r) 
Jo 



x / dt'fi (t') Q (t' - t) 



(20) 



For weak fields, where the exponential can be approxi- 
mated by a first order Taylor series, this expression co- 
incides with the weak field approximation, Eq. (0). In 
analogy to Eq. (fTTj) . we can write the steady state inver- 
sion as 



W s = Wq exp 



dw \Q{uj)\ 2 H (w) 



(21) 



with the spectral filter function H (w) defined in Eq. (|T2j). 
and the Fourier transform f2 (uj) given in Eq. (|10|) . As dis- 
cussed in Section 4, this expression contains a CE phase 
sensitivity only for pulses with phase dependent power 
spectrum. If 72 is large compared to oj(, a , w c and the 
spectral width of the electric field, the steady state in- 
version is given by 



w s = wq exp 

= Wq exp 



4 f°° 
/ duj \Q (oj) 

""72 Jo 

4 f°° 

— / dt \n{t)\ 

72 J-00 



(22) 



which does not depend on the CE phase or the pulse 
shape, but only on the total pulse energy. This result 
can also be obtained by adiabatic elimination of the po- 
larization in Eq. (U). 

In Figs. [5] - QTJ the results of the numerical simula- 
tion (solid curves) and the strong dephasing approxima- 
tion (dashed curves) for the CE phase dependent steady 
state inversion are shown as a function of 72 /lo c . In 
Fig. [8j the phase averaged inversion w s and the mod- 
ulation amplitude A are displayed for interaction with a 
rectangular pulse, assuming a moderate Rabi frequency 
Or = 0.25 uj c . The results of the approximation and the 
numerical simulation agree well, especially for increased 
dephasing, 72 > 0.1 uj c . Fig. Oshows the phase averaged 
inversion w s for excitation by sinc-shaped pulses with the 
same parameters as before. Here, even better agreement 
between numerical and analytical results is found. The 
CE phase dependence of the inversion is negligible in this 
case, |A| < 10~ 7 , since nonlinear phase sensitive effects 
do not play a significant role at those field strengths. 
Figs. [TOl and [TT1 show again the results for excitation by 
rectangular and sinc-shaped pulses, respectively, but now 
for a higher Rabi frequency J7r = lo c . In this case, the 
validity range of the strong dephasing approximation is 
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FIG. 8. Excitation with a rectangular two-cycle pulse: Exact 
numerical result (solid curves) and strong dephasing approx- 
imation (dashed curves) for the (a) average inversion w s and 
(b) modulation amplitude A as a function of the normalized 
dephasing rate 72/^- The Rabi frequency is J1r = 0.25 u; c , 
the transition frequency is oJba = 1-5 ui c . 




FIG. 9. Excitation with a sinc-shaped two-cycle pulse: Ex- 
act numerical result (solid curve) and strong dephasing ap- 
proximation (dashed curve) for the average inversion w s as 
a function of the normalized dephasing rate 72/^- The 
Rabi frequency is Ar = 0.25 u c , the transition frequency is 
io ba = 1.5 un- 



reduced as compared to Or, = 0.25 w c , and good agree- 
ment with the exact result is only obtained for 72 > uj c - 
The inversion dynamics under the influence of phase 
relaxation can be classified into different regimes. In the 
Rabi flopping regime, associated with moderate dephas- 
ing, the inversion dynamics is governed by carrier-wave 
Rabi oscillations. As discussed in Section 5, the steady 
state inversion exhibits a significant phase sensitivity for 
sufficiently short and strong pulses, or for pulses with 
phase dependent power spectra. For higher dephasing 
rates, the near-transparency regime is reached, charac- 
terized by an almost CE phase independent steady state 
inversion w s ~ 0. This regime can be associated with 
dephasing values of roughly 72 /w c « 0.5 ... 5 in Fig. [TU1 
and 72 /u! c ~ 0.2 ... 5 in Fig. [TT] It does not exist if the 
excitation is too weak, see Figs. [8jand[9] In the limit of 
strong dephasing, w s approaches its initial value wq, since 
the evolution of the inversion is suppressed. The freezing 
of the quantum states can be interpreted as the quan- 
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FIG. 10. Excitation with a rectangular two-cycle pulse: Exact 
numerical result (solid curves) and strong dephasing approx- 
imation (dashed curves) for the (a) average inversion w a and 
(b) modulation amplitude A as a function of the normalized 
dephasing rate 72 /w c - The Rabi frequency is Or = uu c , the 
transition frequency is u>ba = 1.5 UJ C . The inset shows A at an 
increased scale. 




FIG. 11. Excitation with a sinc-shaped two-cycle pulse: Exact 
numerical result (solid curves) and strong dephasing approx- 
imation (dashed curve) for the (a) average inversion w B and 
(b) modulation amplitude A as a function of the normalized 
dephasing rate 72 /w c - The Rabi frequency is J1r = uu c , the 
transition frequency is Uba = 1.5 w c . For A, only the numeri- 
cal result is shown, since the approximation is not CE phase 
sensitive for sine pulses. 



turn Zeno effect p2 where the decoherence of the states is 
induced by frequent measurement or, more generally, by 
rapid system-environment interactions. In the transition 
regime between transparency and freezing of the states, 
the phase sensitivity slightly recovers for the rectangular 
pulse, reaching an extremal value of A = —1.8 x 1CP 3 
in the inset of Fig. 1101 In the quantum Zeno regime, 
corresponding to 72/^0 ^ 5 x 10 3 in Figs. HDlandHH the 
modulation amplitude A and also the modulation depth 
6 = A/ (w s — wq) asymptotically approach zero. 



CONCLUSION 



In conclusion, we have studied the carrier-envelope 
phase sensitivity of the steady state inversion w s , re- 
maining in a two-level system after interaction with a 
pulse. We examined the invariance of w s under time re- 
versal of the field. Based on this property, we introduced 
as parameters the phase averaged inversion w s and the 
modulation amplitude A. We discussed the two-level dy- 
namics analytically in the weak field limit, where a phase 
sensitivity of the inversion could only be found for pulses 
with phase dependent power spectrum. Beyond the lin- 
ear regime, the phase sensitivity of w s was studied nu- 
merically for excitation with rectangular and sinc-shaped 
pulses. For the rectangular pulse, it was furthermore 
proven that the modulation depth 8 = A/ (w s — wq) does 
not depend on the pulse length, and an approximate an- 
alytical expression for S was derived applying almost de- 
generate perturbation theory. The influence of phase re- 
laxation on the inversion was investigated based on the 
strong dephasing approximation and numerical simula- 
tions. 



For interaction with few-cycle laser pulses, the two- 
level model predicts a considerable phase dependence 
of w s if the Rabi frequency approaches or even exceeds 
the carrier frequency, see Fig. [7] The necessary field 
strengths can be reached even for unamplified femtosec- 
ond pulses directly out of the laser oscillator by tightly 
focusing the light onto the sampled These model cal- 
culations show that the phase dependent inversion is a 
promising route for the CE phase detection of few-cycle 
laser pulses. One has to be aware of the fact that the ap- 
plicability of the two-level model to an atomic, molecular 
or solid state system in strong fields has its limitations. 
Nevertheless, recent theoretical and experimental work 
shows that basic features of the two-level model are pre- 
served for the strong field excitation with few-cycle laser 
pulsesi 18 ' 19 



The generated inversion can be detected optically, ei- 
ther by an additional probe beam or by the short wave- 
length spectral component of the pulse itself, which is 
delayed against the excitation part of the pulse and can 
for example be chosen at resonance with the two-level 
system. It is also possible to read out the excited elec- 
trons electronically with a field ionization pulse. 8 Exper- 
imentally such a detector could be realized by a gas or 
a metal vapor. For example, in the two-level approxima- 
tion, we obtain with Or, = 1.5w c a modulation amplitude 
of A = —0.34 for interaction of two-cycle Ti:sapphirc 
laser pulses with a transition around 589 nm, correspond- 
ing to the famous D-lines in sodium. Another option 
might be the use of artificial quantum systems like quan- 
tum dots. 
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APPENDIX A: INVERSION INVARIANCE UNDER TIME 
REVERSAL OF THE FIELD 



out the last line of Eq. (|A3[) . we see that the solution 
is only invariant under time reversal of an arbitrary field 
Q (t) if 



M_ ni M_„ 2 . . . M_ ne u = M ne ... M„ 2 M ni u (A4) 



We now want to show that, for 71 = and vanishing 
initial dipole moment, si = s 2 — for t — > — 00, the 
remaining inversion after passage of the pulse is invariant 
under time reversal of the field ft (t) . This should not be 

confused with the time-reversed dynamics of the system. for any subset of indices n with —N < m < n 2 < . . . < 
The Bloch equations in the rotating frame, Eq. (gj, can ni<N. For an even number i, we obtain 
in matrix form be written as 



u = 2f7Afu- 72 (w 1 ,M 2 ,0) ] 



(Al) 



with 



M(t) 



M ne ... M n2 M ni u = w {-if 2 

1/2 

x n 

COS [wfca {tn 2m -i '»2m)] 



sin (uj ba t) 

cos (ujb a t) I (A2) 

■sin(w fca t) -cas(u ba t) 



x(o,o,ir 



(A5) 



and u = (ui,U2,u 3 ) . We assume that the pulse extends 
from —to to to. As initial condition at t — —to, we choose 
u = (0,0,w ) T - First we assume 72 = 0. The formal M„ 
solution of Eq. (|A1[) is then given by 



u (t ) = Texp 



If £ is odd, we obtain 



2fl (t) M (t) dt u 



.M n2 M ni u = w (-lf~ 1}/2 

[t-l)/2 

X Yl C0S [^ba (tn 2m -i ~ tnam)] 
m— 1 

x (sin(ui ba te) , cos(w bQ t^) , 0) T .(A6) 



lim TT exp[20(t_ n )M_ n t /iV]uo(A3) 

I— >-oo 

n=-N 
N 

lim TT [I + 20 (t- n ) M- n to/N] uo 



Thus, for odd i, Eq. (TA4|) is only fulfilled for the third 
vector component, and thus only S3 (to) is invariant under 
time reversal of the field. 
n= _ N The case with 72 ^ can be handled analogously 

by using the transformation u\ = u\ exp (— 72O, U2 = 
with the time ordering operator T, the unity matrix 7, U2 exp (— 72t). Eq. (|A4[) must now be fulfilled for the 
t n = tgn/N and M n = M (t n ). By formally multiplying matrix 

I 



M(t) 



- sin (uj ba t) exp (-72*) - cos (u ba t) exp (-72*) 

I 

For an even number £, we obtain 



sin (ubat) exp(7 2 t) 
cos (w fca t) exp (7 2 t) 




M„ 



M n2 M ni u = wo (-1) 



e/2 



x JJ {exp [72 (t„ 2m _ 1 - t n2m )] 

771=1 

x cos [uj ba (t„ 2m _ 1 - t„ 2m )] } 
x(0,0,l) T . (A8) 



(A7) 



If £ is odd, we obtain 

M nt ... M n2 M ni u = wo (-lf- 1)/2 
(e-i)/2 

X IT { eXp t 72 (*™=m-i ^ *nam)] 



X COS [w ba (t„ 2m _ 1 - t„ 2m )] } 

x exp (72^) 

x (sin (oj ba ti) , cos {uj ba t e ) , 0) T .(A9) 



Again, Eq. (IA4[) is only fulfilled for the third vector com- 
ponent, and thus only S3 (to) is invariant under time re- 
versal of the field. 
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APPENDIX B: MODULATION DEPTH FOR THE 
RECTANGULAR PULSE 

In the following, we want to show that for interac- 
tion with a rectangular pulse shape, the modulation 
depth S, defined in Eq. (IT6l) . does not depend on 
the pulse duration T. Assuming no dissipation, we 
can describe the state of the two-level system by the 
ket (t)} = c a (t)\a) + Cb (t) \b), with the state vec- 
tor c (t) = (c a (t) , Cb {t)) T . The initial condition p a b — 
(si + 1S2) /2 = c*q, = 0, see Section 3, is here imple- 
mented by choosing the initial state vector co=(l,0) , 
corresponding to wq = — 1. The evolution of the system 
can be described by the Schrodinger equation. For an 
exciting field with period T , we can bring the solution 
into the form 

c„ = U n co, 

where U is the unitary transformation matrix, and c„ is 
the state vector after n periods^ We diagonalize U by 
means of a unitary matrix T, c' n = Tc„, U' = TUT^ 1 , 
and obtain for the jth component of cl 



population ppp (t) = [w (t) + 1] /2 in a two-level system 



exp (-inWjTo/h) c/ 0j . 



Here, Wj are the Floquet energies, which do not depend 
on the CE phase 4>ce of the exciting field. By back trans- 
formation, we obtain the occupation probability of the 
upper level 



\c b (nT )\' 



[cxp(-inW b T /h)-cxp(-inW a To/h)]Ta.. 
T aa T bb ~T ab T ba 



(Bl) 



and the inversion is given by w (tiTq) = 2 |c& (nTb)| — 1. 
The matrix elements Ty depend on </>ce, but not on the 
time uTq. Thus, the time dependence cancels out in the 
expression for the modulation depth, Eq. p^|) . We note 
that for the initial condition si = S2 = 0, the modulation 
depth is invariant with respect to wq, since the choice of 
wq just leads to a scaling of the Bloch vector solution s (t) 
in Eq. ((T|). Thus, the proof is also valid for wq 7^ — 1. 



APPENDIX C: ALMOST DEGENERATE 
PERTURBATION THEORY 

Following the formalism in Ref. 28, we can derive a 
perturbative result for the phase dependent inversion. In 
this Appendix, the references to equations refer to Ref. 
28, unless otherwise stated, and we adopt the nomencla- 
ture used there, with wo = <^ba, w = w c , and A = — Or./2. 
As initial condition, we assume wq = —1 as in Ap- 
pendix B (also see remark at the end of Appendix B), 
corresponding to p aa (t ) = 1, ppp (to) = p a p (to) = 
Ppa(to) — in Ref. 28. From Eq. (1.33), we get af- 
ter a few straightforward manipulations the upper-state 



P/3/30) 




m 

2 

^(a,m| X „)e i " mt ° 

m 

a,m\ X+ )e [ " mt ° 



(CI) 



Here, \a,m), \/3,m) are the pseudostates, and |x+), \x~) 
are the Floquet modes of the two-level system. The fre- 
quency q is given in Eq. (1.71). With Eqs. (1.66), (1.73), 
(1.74), (1.76), we obtain in first-order approximation 

\ X+ )=N[ C os(6)(\a,a)-r,\l3,l)) 

+ sin(0)(|/?,-l)+r?k-2))], 

\ X -) = N[-sm(9)(\a,0)-r)\p,l)) (C2) 
+ cos (6) (|/3,-l) + r)\a,-2))], 

with the normalized field rj = A/ (uj + ujq). The normal- 
ization constant N is given in Eq. (1.80b), and sin (6), 
cos (8) are defined in Eq. (1.72). Furthermore, by in- 
serting Eq. (1.77) in Eq. (1.73), we obtain within the 
first-order approximation 

ejj = —haj /2 — rfh (cj + w) , 

efefc = frojo/2 — hu> + rj 2 h (uj + lj) , (C3) 

£jk = tkj = i)h (wo + w) - rfh (uj + uj) . 

The case </>ce = corresponds to a cos pulse extending 
from to = —T/2 to t = T/2, where T is equal to or an 
integer multiple of the carrier period T = 2it/uj. We 
obtain for the remaining population after passage of the 
pulse 

ppp (T/2) = 2N 4 [sm(9) - V cos (6)f 

x [cos(0) +rjsm(6)f [1 - cos(2gT)] (C4) 

The case <^ce = — tt/2 corresponds to a cos pulse ex- 
tending from t = -T/2 - T /4 to t = T/2 - T /4. We 
obtain 

ppp (T/2 - T / A) = 2N 4 [sin (6) + r, cos (B)f 



x [cos (9) — 77 sin (0)]" 
x [1 - cos (2qT)] . 



(C5) 



12 



With Eq. (|T6|) in our paper, we get the modulation depth 



S = 



2rj sin (0) cos (9) sin (9) 2 - cos (6>) 2 (l - ?/ 2 ) 



sin(0) 2 cos((9) 2 (1 



6?7 2 + t] 4 ) + rj 2 



(C6) 

We now insert Eq. (1.72), making use of the relations 
sin 2 (<9) - cos 2 (9) = -A/q, sin (0) cos (0) = e jk / (2qH). 
The frequencies and A are given in Eq. (1.71). We mul- 
tiply the enumerator and denominator of our expression 
for S by g 2 /?y 2 , and subsequently keep only the terms up 
to second order in r\. Making the transition to the nomen- 
clature used in our paper, with ujq — > Wb a , uj — > uj c , and 
7] — > — Qr/ [2 (uj ba +ui c )], we finally obtain the perturba- 
tive result given in Eq. (|17[) in our paper. 
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